Electron-Magnon Scattering in Anomalous Hall Effect 
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We study the role played by electron-magnon scattering in the anomalous Hall effect. We find 
that it has important contributions distinct from other scattering processes like impurities scattering 
and phonon scattering. As a demonstration, we calculate the Hall conductivity for a two dimen- 
sional Dirac model. The result indicates that as system control parameter varies, the competition 
between magnon scattering and other types of scattering changes the Hall conductivity drastically. 
In particular, the side jump contribution could acquire a strong temperature dependence. 
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The anomalous Hall effect (AHE), in which a trans- 
verse voltage is induced by a longitudinal current flow 
in ferromagnetic materials, is one of the most intriguing 
effects in physics. While it has been widely used experi- 
mentally as a standard technique for the characterization 
of ferromagnets, the theoretical study of the AHE proves 
to be complicated and is a subject full of controversial 
issues and conflicting results [lj. In recent years, an 
important connection has been established between the 
AHE and the Berry phase of Bloch electrons @, H, 0, 
This triggers revived interest in this subject and is fol- 
lowed by extensive researches both theoretically and ex- 
perimentally It is now generally accepted that the 
AHE is due to spin-orbit coupling in ferromagnets, and 
apart from an intrinsic contribution which is scattering 
independent, there are also important extrinsic contri- 
butions to the AHE. Especially, there is a peculiar side 
jump contribution that arises from scattering, but does 
not depend on the scattering strength. 

Up to now, most of the theoretical studies of the AHE 
only take into account the impurity scattering, despite 
that many experimental measurements are performed at 
finite temperatures hence other scattering processes like 
phonon scattering and magnon scattering should also be 
relevant. It has been shown that the phonon scattering 
produces similar contributions as impurity scattering 0] , 
which seems to imply that the side jump contribution 
should only weakly depend on temperature. However, 
recent experiment by Tian et al. shows that the side 
jump does have strong temperature dependence Q. 

In this paper, we show that the magnon scattering, 
which has been largely overlooked so far, has distinct 
contributions to the AHE from both impurity scattering 
and phonon scattering. The underlying reason for this 
difference is that the magnon scattering involves spin-flip, 
whereas the impurity scattering and the phonon scatter- 
ing are both spin independent. The extrinsic contribu- 
tion to the AHE turns out to depend sensitively on the 
type of spin structure of the scattering process, which 
has also been noticed in the recent study of the mag- 
netic impurity scattering 0, E3]. Therefore, the anoma- 



lous Hall conductivity would change drastically as two or 
more types of scattering compete. In particular, the side 
jump contribution will depend on the relative scattering 
strength of different types of scattering processes, hence 
should have strong temperature dependence. 

We begin with a description of the electron-magnon 
scattering process in ferromagnetic systems. The ex- 
change coupling between conduction electron and local 
magnetic order parameter can be written as 



flint = - J J dT i & ( r ) ■ S ( r ) 
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where J is the exchange coupling constant, <r is the vector 
of Pauli matrices for conduction electron spin, S is the 
local spin, a± = & x ± ia y , S± = S x ± iS y , and hat means 
the quantity is a matrix in spin space. The last term 
above describes the exchange splitting which should be 
included in the non-interacting part of the Hamiltonian, 
whereas the first two terms describe the electron-magnon 
scattering. Using Holstein-Primakoff representation, for 
temperatures below Curie point when the number of local 
spin flips is much smaller than the total spin, we can write 
the interaction Hamiltonian as 
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where V is the system volume, c 1 (c) and (a) are 
the electron and magnon creation (annihilation) oper- 
ators respectively. From Eq.([2]), it is clear that the 
magnon scattering flips spin, unlike impurity scattering 
and phonon scattering which are spin independent. 

When the energy of magnons involved in the scatter- 
ing is much lower than the Fermi energy and the electron 
energy spectrum (and density of states) varies smoothly 
around Fermi surface, we can approximate the scattering 
process as quasi-elastic. In typical ferromagnetic materi- 
als like Fe or Co, the magnon has a large effective mass 
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about 1CT 29 ~ 10 _28 kg [Tl(. The quasi-elastic treat- 
ment will be a good approximation if the electron effec- 
tive mass at Fermi level is much smaller than that of 
the magnon. In this case, the electron sees an effective 
scattering potential 



(3) 



where V£(q) = —Jy/Sn m (q)/2 is the orbital part of the 
scattering potential, n m (q) is the distribution function of 
magnon, and the spin part <r± is off-diagonal representing 
spin-flip processes. 

To demonstrate that the magnon scattering gives dis- 
tinct contributions to the AHE, we calculate the Hall 
conductivity of the two dimensional (2D) Dirac model. 
The AHE from impurity scattering in this model has 
been studied previously by Sinitsyn et al. [13] . We choose 
this model not only because of its simplicity to demon- 
strate our ideas, but also because it describes low en- 
ergy physics of interesting systems such as graphene [3] , 
kagome lattice 14| and surface states of topological insu- 
lator [15[ (though it should be noticed that for graphene 
and kagome lattice the 'spin' refers to the sublattice de- 
grees of freedom rather than electron spin as discussed 
here). Therefore the results presented here will also be 
important in understanding transport properties of these 
systems. 

The AHE occurs in 2D Dirac model when a suitable 
symmetry breaking mechanism is introduced. Our model 
Hamiltonian reads (we set ft = 1 in the following) 



H = v(k x a x 
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where the last term is the symmetry breaking term which 
introduces a gap of 2A. We assume that this model de- 
scribes low energy physics near the Fermi surface of cer- 
tain 2D ferromagnetic system, hence the term Aa z rep- 
resents the exchange splitting which corresponds to the 
last term in Eq.JT]). 

The eigenstates of the system are given by "0fc V ( r ) = 
(l/VV)e lkr |w^,' v ) with energy eigenvalues e c,v (fc) = 
±\/ (vk) 2 + A 2 , where c and v stand for conduction and 
valence band respectively, and \u^ v ) is the spin part of 
the eigenstate which can be written as 



sin ; 
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A/^(vk) 2 + A 2 and 



y /k x . Due to spin-orbit coupling, the eigen- 



with 9 and </> being the spherical angles of the vec 
tor (vk x , vk y , A) such that cos0 
tan0 = k 
spinors are fe-dependent. 

We evaluate the Hall conductivity by using the Kubo- 
Streda formalism [H E3. In this approach, the Hall 
conductivity can be separated into two parts, a xy = 



and a^y is a Fermi sea contribution for which we only 
need to retain the scattering-free component in the weak 
disorder limit 12]. All the important scattering effects 
are contained in ctL,, which takes the form 
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where G R and G A are the retarded and advanced Green's 
functions respectively, v x<y = va x . y are the velocity op- 
erators, £f is the Fermi energy, the trace is taken over 
both momentum and spin spaces and the bracket means 
statistical average over disorder configurations. In the 
following calculation, we take the Fermi energy to be in 
the conduction band. 

The intrinsic contribution results from the Berry cur- 
vatures of spin-orbit coupled bands and is independent of 
scattering. In Kubo-Streda formalism, the intrinsic con- 
tribution comes from the scattering-free components of 
It can be decomposed into two parts 

- ,, - ".,-„ , where (Txy^ is the contribution 
from all the completely occupied valence bands below 
the Fermi surface and cr^ is the contribution from the 
partially filled conduction band where the Fermi surface 
lies in. The contribution from completely filled bands 
a X y has a topologically quantized value Ne 2 /{2ir) with 
iV" being an integer known as the first Chern number. 
The calculation of N goes beyond any low energy effec- 
tive model since it involves the entire Fermi sea. On 
the contrary, the contribution crj^ from the partially 
filled conduction band can be regarded as a Fermi surface 
property 18|. For Dirac model, we have 
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it ',.',„ where a\ y is a Fermi surface contribution, 
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where dp is the spherical angle 9 evaluated on the Fermi 
surface when k = kp- 

The extrinsic contribution consists of the side jump 
and the skew scattering. In the semiclassical picture, 
the side jump arises from the coordinate shift of a wave- 
packet during the scattering process, and the skew scat- 
tering appears due to the asymmetry of scattering rate 
for higher order scattering processes [l[. It has been 
clarified recently that the skew scattering contribution 
as defined in the semiclassical picture actually contains 
two different parts: a conventional skew scattering part 
with n _1 dependence and an intrinsic skew scattering 
part with n° dependence with n being the disorder den- 
sity [TH. According to its parametric dependence, the in- 
trinsic skew scattering can be included as part of the side 
jump. In Kubo-Streda formalism, the various contribu- 
tions listed above have been identified with different sets 
of Feynman diagrams in the self-consistent non-crossing 
approximation [121 ]. 

Let's first consider a clean system with only magnon 
scattering. In the quasi-elastic approximation, the dis- 
order lines in Feynman diagrams do not carry energy 
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arguments. Since the population of magnon bath is con- 
served in steady state, each disorder line must have a pair 
of <7_|_ and <7_ at the two ends, corresponding to magnon 
emission and absorption processes. Therefore the conven- 
tional skew scattering which involves third order scatter- 
ing events must vanish. Furthermore, due to the angular 
average at velocity vertices, the intrinsic skew scattering 
also vanishes, leaving only the side jump contribution, 

< = ^ cos0 f, (8) 
which cancels with the part of intrinsic contribution that 



depends on Fermi energy, such that the final result be- 
comes a constant value e 2 /(47r) and is the same as the 
intrinsic contribution for a completely filled conduction 
band, i.e. in the limit Of — > tt/2 |20| . 



Next we shall include the impurity scattering as well 
and investigate the competition between magnon scatter- 
ing and impurity scattering in the AHE. For simplicity, 
we consider the short range impurity as been studied in 
Ref. The result of the total Hall conductivity is 

(including only a^y^ for intrinsic contribution) 
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where C = r m 1 / T i 1 ^ s the ratio of magnon scattering rate 
to impurity scattering rate with 
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and j] = t s/c 2 /tj 2 represents the conventional skew scat- 
tering contribution from impurities. For standard white 
noise impurity model, the third order correlation of scat- 
tering potential is zero, so r\ vanishes identically. It will 
be nonzero if a non-Gaussian part V\ of the impurity 



potential is included [l2j, then T sfe 2 = riiE F V± / (4irv 4 ) . 



Observe that when the magnon scattering is dominant 
over impurity scattering, i.e. for very large £, we recover 
the result of Eqs. (|7l8p with a constant value e 2 /(47r). In 
the opposite limit, when impurity scattering dominates, 
C — > 0, we retain the previous result by Sinitsyn et al. [12J. 
From our result Eq.©, it is clear that as ( varies, which 
results from the competition between different scattering 
mechanisms, the extrinsic contributions to the Hall con- 
ductivity varies drastically and can have a sign change. 

It is of fundamental importance to experimentally sep- 
arate contributions from different mechanisms, especially 
the intrinsic contribution Q. Let us collect the terms of 
a xy that are of order n , denoted as a® y . These include 
the intrinsic contribution and the side jump (including 
intrinsic skew scattering). For soft magnon modes, the 
ratio C does not sensitively depend on the Fermi energy. 
In Fig|TJ we plot a® y as a function of Fermi energy for 
different values of £. As ( increases, the curve of Hall con- 
ductivity is shifting upward from the impurity scattering 
dominated situation and approaching the limiting value 
e 2 /47r for the magnon scattering dominated case. This 
competition behavior is more clearly observed in Figj^l 
where a® y is plotted at fixed Fermi level as a function 
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FIG. 1: (color online). a" y plotted as a function of the Fermi 
energy ef for fixed values of £. a° y is measured in units of 
e 2 /(47r), and ef is measured in units of A which is half of the 
gap size. 



of £. As C increases, a® y increases monotonically. For 
impurity scattering dominated case, a xy takes negative 
value for Fermi energies below ef ~ 7. 3 A. Hence in this 
energy range, there is a sign change of a® as C increases, 
i.e. when magnon scattering gradually takes dominant 
place. 

As observed from this model calculation, the magnon 
scattering indeed plays a quite different role as compared 
with impurity scattering. This difference comes from 
their different structures in spin space. Magnon scatter- 
ing flips spin hence its the matrix element is off-diagonal 
in spin space while both impurity scattering and phonon 
scattering is proportional to the identity in spin space. 
As a result, the extrinsic contributions for each type of 
scattering involve different combinations of structure fac- 
tors (such as cos Of and sin Op in our example), which 
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As a result, the side jump contribution can have strong 
temperature dependence. 

The authors thank A. H. MacDonald and Y. Yao 
for valuable discussions. This work is supported 
by DOE (DE-FG03-02ER45985), NSF(DMR0906025), 
Welch Foundation (F-1255), and NSFC (10740420252). 



FIG. 2: (color online). a® y versus the ratio of scattering rates 
(" for fixed values of Fermi energy ef- cr^y i s measured in 
units of e 2 /(47r). The plot shows the crossover from impurity 
scattering dominated regime to magnon scattering dominated 
regime as £ increases. 



leads to the competition behavior. Magnetic impurity 
scattering has yet another spin structure different from 
the above three, as being proportional to a z if the aver- 
age magnetization is along z-direction. Previous studies 
indeed show that the magnetic impurity scattering be- 
haves differently from the normal impurity scattering for 
the AHE [^, [l(|. The above analysis suggests that we 
could classify various scattering processes according to 
their structures in spin space which is the major fac- 
tor that determines their contributions to the AHE 21]. 
The competition between different classes could change 
the Hall conductivity dramatically as system control pa- 
rameter such as temperature varies. 

Finally, we point out that the above discussion is not 
limited to ferromagnetic systems with electron spin de- 
grees of freedom. Any quasi-particle index which has 
two degrees of freedom can be generally referred to as 
'spin' (or pseudospin). Anomalous Hall effect will arise if 
the system has 'spin'-orbit coupling as well as 'spin' split- 
ting. For example, in a bipartite lattice such as graphene, 
the sublattice degree of freedom can be treated as pseu- 
dospin. Anomalous Hall transport occurs in graphene 
when there is sublattice symmetry breaking in the sys- 
tem (22| . As another example, for bilayer systems, it is 
the layer index that plays the role of pseudospin and the 
pseudospin splitting can be realized by imposing a bias 
between the two layers. In general, our result indicates 
that a careful analysis of various scattering processes ac- 
cording to its pseudospin structure is indispensable in the 
study of AHE for these systems. 

In summary, we have shown that the electron-magnon 
scattering plays an important role in the anomalous Hall 
effect. The competition between magnon scattering and 
other scatterings changes the anomalous Hall conductiv- 
ity drastically as system control parameters are varied. 
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